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j^ I Abstract 

_L ' The design of vesicles which become unstable at an easily tuned value of pH is of 

Q . great interest for targeted drug delivery. We present a microscopic theory for two forms 

o : 

of such vesicles. A model of lipids introduced by us previously is applied to a system of 

^j>. I ionizable, anionic lipid, and permanently charged, cationic lipid. We calculate the pH 

\^ ' at which the lamellar phase becomes unstable with respect to an inverted hexagonal 

(N' 

On ■ one, a value which depends continuously on the system composition. Identifying this 

O ; 

^^ . instability with that displayed by unilamellar vesicles undergoing fusion, we obtain very 

"^ , good agreement with the recent experimental data of Hafez et al., Biophys. J. 2000 
79: 1438-1446, on the pH at which fusion occurs vs. vesicle composition. We explicate 

pj ' the mechanism in terms of the role of the counter ions. This understanding suggests 

o' 

O ■ that a system of a neutral, non lamellar forming lipid stabilized by an anionic lipid 



X 



would serve equally well for preparing tunable, pH sensitive vesicles. Our calculations 
confirm this. Further, we show that both forms of vesicle have the desirable feature of 
exhibiting a regime in which the pH at instability is a rapidly varying function of the 
vesicle composition. 



I. Introduction 

The creation of liposomes which become unstable in response to changes in their environment 
has been the object of longstanding interest in connection with applications to drug delivery 
(Yatvin et al., 1980). One particularly interesting environmental cue is the relatively low 
pH which is found in tumor tissue (Tannock and Rotin, 1989) and in endosomes (Tycko and 
Maxfield, 1982). In the latter case, the rapid acidification which occurs in the endocytic 
vesicle would bring about the instability of the liposome resulting either in the release of its 
contents within the endosome itself or in liposome-mediated destabilization of the endocytic 
vesicle with consequent release of the liposome's contents to the cytoplasm (Straubinger, 
1993). 

There exist various strategies for producing pH-sensitive liposomes (Thomas and Tirrell, 
1992; Torchilin et al., 1993; Straubinger, 1993; Chu and Szoka, 1994; Sorgi and Huang, 1996). 
One method is to combine a lipid which does not form bilayers under physiological conditions 
with an ionizable anionic amphiphile or lipid. The latter, when sufficiently charged, stabilizes 
a bilayer of the combined system. The mechanism of this stabilization, as we argue below, 
is the attraction of counter ions and their associated waters of hydration to the vicinity of 
the headgroups which effectively increases their size. As the pH is reduced, so is the fraction 
of anionic amphiphiles which are ionized. Therefore there are fewer counter ions near the 
headgroups to stabilize them. Thus the reduction in pH eventually triggers an instability of 
the vesicle; the lipids revert to their more stable phase, usually an inverted hexagonal (Hji) 
one. Most commonly, the non lamellar-forming lipid is a phosphatidylethanolamine (PE), 
such as dioleoylphosphatidylethanolamine (DOPE) (Cullis and de Kruijff, 1978). 

One problem with this strategy is that the pH at which the instability occurs is deter- 
mined by the pK of the single ionizable component, and is therefore not easily tuned. Discrete 
tuning can be obtained by utilizing different ionizable components (Collins et al. 1989). An 
alternative method which results in a vesicle exhibiting an instability at a value of pH which 
can be tuned continuously has recently been demonstrated (Hafez et al., 2000). They utilized 



a vesicle with an ionizable anionic lipid, cholesteryl liemisuccinate (CHEMS), and a perma- 
nently charged cationic lipid N,N-dioleoyl-N,N-diniethylammonium chloride (DODAC). The 
pH at which the vesicle becomes unstable is a monotonically increasing function of the DO- 
DAC concentration, and is therefore easily, and continuously, tuned. We understand this 
result as follows. Vesicles consisting only of CHEMS are unstable with respect to formation 
of an Hu phase at pH less than 4.2 (Hafez and Cullis, 2000). This implies that a sufficient 
number of CHEMS must be ionized to stabilize such a vesicle. These ionized headgroups 
attract counter ions to their vicinity. It is these counter ions, enlarged by their waters of 
hydration, that stabilize a system which would otherwise tend to an Hu phase. The addition 
of cationic DODAC to such a vesicle at any pH, causes a decrease in the number of counter 
ions in the vicinity of the headgroups. This decrease tends to destabilize a previously stable 
vesicle. To restore the number of counter ions near the headgroups, and the vesicle's stabil- 
ity, more counter ions must be attracted to the headgroups, which can be done by ionizing 
more of the CHEMS; i.e. by increasing the pH. Hence the value of the pH at the instability 
is an increasing and continuous function of the concentration of DODAC. 

In this paper, we apply our model to the system of mixed, ionizable, anionic lipid and fully 
charged, cationic lipid and solvent such as that examined by Hafez and Cullis. We identify 
the pH at which bilayer vesicles become unstable as the pH at which the lamellar phase 
becomes unstable to the inverted hexagonal phase, a reasonable assumption supported by 
much experimental data (Hope at al., 1983; Ellens at al., 1986). With this identification, we 
indeed find that the pH at which an instability occurs is a monotonically increasing function 
of the cationic lipid concentration, or equivalently, a monotonically decreasing function of 
the anionic lipid concentration. Our results fit the experimental data very well. 

The above reasoning indicates that tunable, pH sensitive liposomes should also be formed 
from a mixture of a neutral lipid which favors a non-lamellar phase, such as phosphatidyl- 
ethanolamine (PE), and an anionic lipid which stabilizes the liposome by attracting coun- 
terions to it. Such stabilization is well known utilizing various anions, such as palmitoylho- 
mocysteine (Yatvin et al., 1980; Connor and Huang, 1985), oleic acid (Straubinger et al.. 



1985; Wang and Huang, 1987), or CHEMS (Straubinger, 1993; Ellens et al., 1984). As the 
concentration of the stabilizing counter ions clearly depends on both the concentration of 
the anionic lipid and the pH, the value of the latter at which the vesicle becomes unstable 
will be tunable, depending continuously on the concentration of anionic lipid. To test this 
hypothesis, we apply our model to a system of mixed, ionizable, anionic lipid and neutral 
lipid. We again find that the pH at which an instability occurs is a monotonically decreas- 
ing function of the anionic lipid concentration. Our results are in accord with experiments 
on the pH sensitivity of vesicles composed of l-Palmitoyl-2-oleoyl-phosphatidyletanolamine 
(POPE) stabilized with tocopherol hemisuccinate (Jizomoto et al., 1994). Lastly our results 
show that both forms of vesicle have the desirable property of exhibiting a regime in which 
the pH at the instability is very sensitive to the concentration of the anionic lipid. This 
will always be the case whenever a minimum amount of one lipid is required to stabilize the 
formation of vesicles by the mixture. 

In the following section we briefiy review our model of charged lipids (Li and Schick, 
2000a) and of lipid mixtures (Li and Schick, 2000b; Li and Schick 2000c). We then present 
our results. First we show the phase diagram of a single, ionizable, anionic lipid, solvent, 
and counter ions, a diagram which shows the transition from a lamellar phase (La) to an 
inverted hexagonal one {Hjj) as a function of pH. We then consider the mixed system of 
ionizable anionic lipid and completely charged cationic lipid, and show the phase behavior. 
There is a transition between La and Hjj phases which occurs at a value of the pH which 
is a function of the concentration of the ionizable anionic lipid. We also show the spatial 
distribution of the various mass and charge densities in the coexisting phases. Lastly we 
consider the system of ionizable anionic and neutral lipids, and present its phase behavior. 
It also shows a transition between La and Hn phases which occurs at a value of pH that 
depends on the concentration of the ionizable lipid. 



II. The Model and its Self Consistent Field Solution 

We consider a system of volume V consisting of anionic lipids, cationic lipids, counter ions 
and solvent whose densities are controlled by the fugacities Zi, Z2, Zc, and Zg respectively (Li 
and Schick, 2000a). By taking the charge of the cationic lipid to zero, we can also describe 
a mixture of anionic and neutral lipids. The counter ions are positively charged. Because 
of overall charge neutrality, the average amount of charge on the anionic lipids is related to 
the density of cationic lipids and the density of counterions. Hence we use the fugacity Zc to 
control the charge on the anionic headgroups, which is equivalent to controlling the pH. 

With the exception of their Coulombic properties, the two lipids are modelled identically. 
They each consist of headgroups of volume Uh, and two equal-length, completely flexible tails 
each consisting of A^ segments of volume z/^. Each lipid tail is characterized by a radius of 
gyration Rg = (A^a^/6)^/^, with a the statistical segment length. The counter ions are 
characterized by their charge, +e, and their volume, z/c, while the neutral solvent particles 
are characterized by their volume Ug. 

There are eight local densities which specify the state of the system. We measure them all 
with respect to the convenient density u^^. They are the number density of the headgroups 
of the anionic lipids, i^^^^h ('^)' ^^^ '^^ ^^^ cationic lipids, i^h^^^ (r); the number density of 
the tail segments of each lipid, i^h^^t (i") and i^^^^t (i")5 the number density of the solvent, 
z/;^^$s(r) and of the counter ions z/^^$c(r); and the local charge density of the headgroup of 
the anionic lipid, ev^^Pf^ (r) and of the cationic lipid, ev^^P^^ (r). The local charge density 
of the positive counter ions is simply ez/^^$c(r). Note that all functions $(r) and P(r) are 
defined to be dimensionless. 

The interactions among these densities are of two kinds. First there is a repulsive, 
contact interaction between headgroups and tail segments, and also between solvent and tail 
segments. The strength of this interaction is kTuhX, where k is Boltzmann's constant and T 
the absolute temperature. Second there is the Coulomb interaction between all charges. The 



energy per unit volume of the system, expressed in the natural units A;T/z//j, can be written 
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is a dimensionless measure of the strength of the Coulomb interaction, and e is the dielectric 
constant of the solvent. In addition to these interactions, we impose a local incompressibility 
constraint on the system which models the hard core interactions between all particles. Upon 
defining the volume ratios 7^ = Ug/uh, 7c = vd^h-, and 7^ = 2Nut/uh, the incompressibihty 
constraint that the sum of the volume fractions of all components must be unity everywhere 
takes the form 
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As shown earlier (Li and Schick, 2000a), the partition function of the system can be 
written in the form in which the eight fluctuating densities, instead of interacting directly 
with one another, interact indirectly via eight fluctuating fields, here denoted W^^ , W} , 
Ul , with L = 1,2, and Wg, U^- Self consistent field theory results when the fluctuating 
fields and densities are approximated by those values which minimize the free energy, fi, of 
the system in the presence of these fields. The free energy to be minimized has the form 
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Here Qc[f^c] is the partition function of a single counter ion of unit positive charge in an 
external potential Uc, 

Qc[Uc] = |rfReexp[-f/,(R,)], (5) 

Qs[W^s] is the partition function of a single solvent molecule in an external field Wg, 

Qs[Ws] = J dRsexp[-Ws{Rs)], (6) 

and QiiW^h , Wf , f/^ ], given below, is the partition function of a single lipid of type L 
in external fields W^^ , W^ , and U^ . Note that a Lagrange multiplier S(r) has been 
introduced to enforce the incompressibility constraint of Eq. |^. The functions Wf^ i ^h 
etc. which extremize this free energy will be denoted by their corresponding lower case 
letters wf^ and (f)\^ etc. 

It is not difficult to see from the form of the free energy Q and that of the energy E of 
Eq. |l] that the fields acting on the different heads and extremizing the free energy are equal, 
Wfi^ (r) = wf^ (r) = u'/j(r), that the fields acting on the different tails and extremizing the free 
energy are equal wl (r) = Wf (r) = Wt{r), and that the fields acting on all charge densities 
and extremizing the free energy are related, uj^ (r) = u}^ (r) = Uc{r) — 'jc^{r) = u{r). Thus 
there are only five independent functions, Wh{r), Wt{r), Ws{t), ■u(r) and ^(r), and these are 
obtained from the five equations 



w,{r) ^ 2xiV^0j^)(r)+^(r), (7) 

L=l 

wt{v) ^ 2xiV[0,(r) + $:0l^)(r)]+7ie(r), (8) 



^,(r) = 2xiV5:0r^(r)+7.e(r), (9) 

L=l 

^ ^ - An] Rl |r - r'| ' ^ ' 



%0s(r) + 7c0c(r) + El^i^'^W +7i#^(r)]. (11^ 



Because the field ^ can be easily eliminated, one deals essentially with four equations. The 
eight densities are all functionals of the above fields except ^ and, therefore, close the cycle 



of self-consistent equations: 
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Note that one of the self-consistent equations, Eq. |10|, is the non-linear Poisson-Boltzmann 
equation, and u{r) is the electric potential. 

With the aid of the above equations, the self consistent, or mean field, free energy, ^2^/, 
which is the free energy function of Eq. | evaluated at the self consistent field values of the 
densities and fields can be put in the form 

fcT / ^ \ 

- ^rnf = [Yl ZLQL[Wh,Wt,u] + ZcQc[Uc\ + ZsQs[Ws\ 

^h \l=1 J 
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where we have chosen J ^{r)dr = for convenience. All of the above is a simple extension of 
the procedure in our earlier paper (Li and Schick, 2000a) with one exception: previously we 
assumed the counter ion to have neglible volume and included the interaction between its 
charge and the dipole of the solvent so that it would attract waters of hydration and gain 
an effective volume. Here we simply assign the counter ion a volume, Uc, so that we need 
not utilize the interaction between charges and solvent dipoles. 

The fact that the anionic lipids are ionizable has the consequence that the fields Wh, Wt 
and u acting on that lipid can be replaced (Borukhov et. al. 1998; Li and Schick 2000a) by 
Whiffy ifj, and where 

w^^lffir) = Whir) - ln[l +p(exp{M(r)} - 1)]. (18) 
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The parameter p is related to the pK of the headgroup and can therefore be related to the 
fugacity Zc of the counter ions by means of the condition of charge neutrality 
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In practice, we use this fugacity to control the fractional charge on the anionic lipid and 
therefore the pH. 

In the system in which the cationic lipid is fully charged, the fields Wh, Wt and u acting 
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on it can be replaced by wj^^ff, Wt, with 
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Note that from Eqs. |T2| and |n| it immediately follows that the number density of the 
headgroup and its charge density, in units of e, are identical. 
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as they should be because the cationic lipid is always fully charged. 

In the other system which we consider, the second lipid is neutral so that 
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and 
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There remains only to specify how the partition function of the lipids is calculated. As in 
our earlier study (Li and Schick, 2000a), one defines the end-segment distribution function 
q^^\r,s) which satisfies the equation 
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Prom this function, one obtains the partition functions of the lipids. 
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the head and tail densities 
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and the charge density of the anionic hpid head 
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The average fractional charge, fc, on the anionic lipid headgroup follows 
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from which the pH relative to the pK of the anionic lipid headgroup is obtained: 

pH = pK + log, J Y^)- ^^^) 

To summarize: there are five self-consistent equations to be solved for the five fields Wh{r), 
Wi(r), Ws{r), u{r), and ^(r). They are Eqs. ^to |TI]. The fields depend on the eight densities 



(f)}^ (r), (j)i (r), pIj (r) with L = 1,2, 0s (r), and 0c (r) which depend, in turn, on these fields. 
The densities are given by Eqs. [1^, |1^, ^ or ^, ^7\, |2^, and ^ Once the fields and densities 
are obtained, the free energy follows from Eq. |T^. 

Instead of solving these equations in real space, we do so in Fourier space in such a way as 
to guarantee that our solution has the symmetry of either the lamellar or inverted hexagonal 
phases (Matsen and Schick, 1994). Comparison of the free energies of these phases tells us 
which is the globally stable one. We do this for different temperatures, lipid concentrations 
and pH, and thereby map out the phase diagram. 

III. Results 

We first consider the system of the single anionic lipid in solvent. The architecture of this 
lipid is characterized by ■jt, the ratio of the volume of its tail groups to the volume of its 
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head. In choosing the value of this parameter, we have been guided first and foremost by the 
requirement that our model lipid exist in the inverted hexagonal phase when its head group 
is neutral and the system is hydrated in order that it model correctly the behavior of CHEMS 
(Hafez and Cullis, 2000). The value we have chosen, 7^ = 2.5 does indeed produce a model 
lipid which exists in the inverted hexagonal phase over a large region of phase space when 
it is neutralized, as is seen below. We note in passing that this is not an unreasonable value 
when compared to that which follows from volumetric data on the non-lamellar forming lipid 
DOPE (Rand and Fuller, 1994), 7^ = 2.94, a value which most likely assigns some of the 
volume of the waters of hydration, those most tightly bound, to the headgroup. But in our 
model, some of the volume of water of hydration should be included in the volume of the 
head group because the only interaction it has with water occurs when it has a net charge. 
The model interactions, therefore, neglect those waters attracted via their dipole moment to 
the charges of a physical, neutral, head group. 

The solvent is characterized by its relative volume 7^ = Vs/i^h = 0.1, close to the ratio 
of 0.096 appropriate to water and a phosphatidylethanolamine headgroup (Rand and Fuller, 
1994; Kozlov et al. 1994). The counter ions are modeled as H9O4 , a reasonable choice 
(Bell, 1959), and are therefore characterized by their relative volume 7c = t^c/i^h = 0.4. The 
strength of the Coulomb interaction is, again, given by the parameter j3*, Eq. |^, which can 
be written as (3* = ^/Li, where ^ = e'^/ekT is the Bjerrum length, and Li = Vh/AuB? is 
a length characterizing the architecture of lipid 1. Using the value of Rg found earlier (Li 
and Schick, 2000a) to be appropriate to DOPE and a Bjerrum length of 7 A appropriate 
for water, we obtain (3* = 27 and have used this value. The phase diagram of the system 
is shown in Fig. 1 as a function of the effective temperature T* = {2xN)^^ and the pH 
relative to the pK of the lipid headgroup. We observe the characteristic transition between 
the inverted hexagonal and lamellar phases as the pH is increased (Hope and Cullis, 1980; 
Bezrukov et al., 1998). The fugacity of solvent here is Zg = 3.2. At this value, the system 
is not in the presence of excess water. When completely neutralized, the two phases coexist 
at T* = 0.065 with the Hfj phase containing a volume fraction of solvent 7s0s = 0.063 and 

11 



the lamellar phase a volume fraction 7s0s = 0.087. For a given temperature, the transition 
between phases occurs at a given pH. If this value is not a biologically useful one, as it is 
not for PS which undergoes a phase transition at the very acidic value pH ^ 3 (Hope and 
CuUis, 1980), then a vesicle made from this lipid is not applicable for drug delivery. 

To vary continuously the value of the pH at which the lamellar phase becomes unstable 
for a fixed temperature, one can add an additional lipid. We now consider the case when 
this additional lipid is cationic, and fully charged as in the experiment of Hafez et al. (Hafez 
et al., 2000). The phase diagram we obtain for this system is shown in Fig. 2 as a function 
of the fractional concentration of the ionizable anionic lipid 

e = tf'+;'-^i",„. (32) 

VS. pH — pK, where the pK is that of the anionic lipid headgroup. Our results for the 
phase coexistence are shown in solid lines. We have chosen the solvent activity Zg = 3.425 
and the effective temperature to be T* = 0.079. Under these conditions, the system of 
pure anionic lipid is, when completely neutralized, just at phase coexistence between L^ 
and Hji phases. Therefore our curves, in the limit of very large negative pH, asymptote to 
G = 1. This coexistence is characterized by a volume fraction of solvent 7^0^ = 0.087 in 
the inverted hexagonal phase and 7s0s = 0.116 in the lamellar phase. Because the phase 
boundary curves are flat near = 1.0, our results are not very sensitive to different choices 
of temperature and solvent chemical potential provided, of course, that we remain in the 
same general region of phase behavior. The experimental results of Hafez et al., which show 
the pH at which vesicle fusion occurs (their pH/), are shown by the solid dots. To obtain the 
best fit, we have taken the pK of CHEMS to be 5.5. The actual value of the pK of CHEMS 
in the DODAC/CHEMS system has not been measured. It has been measured in a large 
unilamellar vesicle composed of CHEMS and DOPE and a value of 5.8 obtained (Hafez and 
Cullis, 2000). Thus the value of 5.5 we have used in our fit to the data is not unreasonable. 
We note that our results fit the data rather well except for the very largest value of pH. 
Given the assumptions in the modelling, and the assumption that the architecture of the 
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two lipids is the same, this agreement is gratifying. We also restate the interesting point 
made by Hafez et al. that the preferred phase of the lipid mixture can be inverted hexagonal 
even though both of the lipids in isolation adopt a lamellar organization. They do so, in 
our view, because in isolation their headgroups are sufficiently charged to attract stabilizing 
counter ions which increase the effective size of their headgroups. The transition to the 
inverted hexagonal phase comes about because, as the lipids are mixed, the number of those 
stabilizing counter ions is reduced until eventually the lamellar phase becomes unstable. 

A somewhat different point of view is taken by Hafez et al. They assume that for 
fusion to occur, i.e the instability, the surface charge on the vesicle must be zero, permitting 
close contact. Therefore the proportion of CHEMS which is negatively charged must equal 
the DODAC content of the membrane. Equivalently, the number of counter ions must be 
identically zero. This condition is expressed in the equation 

pH-pK = logio (^0^) ' (33) 

which is plotted in Fig. 2 as the dashed line. 

The spatial distributions of the various components and of the charges is shown in Fig. 
3 for the lamellar phase and Fig. 4 for the inverted hexagonal phase which coexist near 
pH— pK=0. In part (a) of each figure, the volume fractions of all elements are shown: of the 
headgroup of the anionic lipid, (pj^ , of the tails of the anionic lipid, 'jt<Pt > of the headgroup 

(2) (2) 

of the cationic lipid, (f)\ , and of the tails of the cationic lipid, 7(0^ , of the solvent, 7s0s, and 
of the counter ions, 7c0c- In part (b) of each figure, the charge distributions, in units of euj^^, 
are shown: that of the anionic lipid, pj^ , of the cationic lipid, pj^ , of the counter ions, 0c, 
and the total charge distribution. It should be recalled that these distributions are those of a 
lamellar phase, not an isolated lipid bilayer. In the latter, the volume fraction of headgroups 
would fall rapidly in the solvent-rich regions on either side of the bilayer. The points x = 
and X = D correspond to the centers of sequential solvent regions within the lamellae, with 
D = 2.94Rg being the lamellar period. In Fig. 4, x = and x = D correspond to the centers 
of adjacent tubes with D = S.OSRg the lattice constant of the inverted hexagonal phase. As 
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noted previously (Li and Schick, 2000a), a single dielectric constant has, for simplicity, been 
employed for the entire system. Were a different dielectric constant employed in the tail 
region, the distribution of counter ions would change, with less of them in the tail region. 
However their volume fraction in that region is already small due to their non-zero volume 
and the incompressibility constraint. Hence any change in the distribution would probably 
be small. 

We turn now to the system of anionic and neutral lipids. The phase diagram is shown in 
Fig. 5. As in Fig. 2, is the fractional composition of the anionic lipid. The temperature 
T* and solvent fugacity are the same as in Fig. 2. The results for the phase coexistence 
of the anionic, neutral lipid system are shown with solid lines. They are compared to the 
calculated results for the anionic, cationic system shown previously in Fig. 2, and repeated 
here in dashed-dotted lines. We see again in this system of anionic and neutral lipid the 
characteristic transition from Hjj to La phases with increasing pH. Again the pH at which 
the transition occurs is a continuous function of the system's composition, decreasing with an 
increase in the composition of the ionizable anionic lipid. This is in agreement with results 
on vesicles of POPE and CHEMS (Jizomoto et al., 1994). Equally important, we note a 
region in which a very small change in the composition of the system brings about a large 
change in the pH at which the instability of the lamellar phase occurs. We also note that a 
minimum fractional composition of anionic lipid, approximately 0.3, is necessary to stabilize 
the DOPE-like neutral lipid. This is similar to the experimental observation that a minimum 
molar composition of 0.2 CHEMS was necessary to stabilize vesicles of transesterified egg 
phosphatidylethanolamine (Lai et al., 1985). It is easy to see from Fig. 5 that the existence 
of a minimum concentration of stabilizing anionic lipid guarntees a regime in which the pH 
at the instability changes rapidly with concentration. We understand the instability in this 
system in the same manner as for the previous system. The neutral lipid prefers to be in 
an inverted hexagonal phase. The addition of an anionic lipid can stabilize a lamellar phase 
of the mixture if there is enough of it, and if these lipids are sufficiently charged. It does 
so by attracting a sufficient number of counter ions which effectively increase the headgroup 
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of the anionic lipid. When the pH is changed so that fewer counter ions are attracted, the 
lamellar phase is less stable, and eventually becomes unstable with respect to the inverted 
hexagonal phase. Note that the criterion of Hafez et al. that the density of counter ions be 
zero at the transition yields no useful information for this system. 

IV. Discussion 

We have applied our model of lipids, introduced previously (Li and Schick, 2000a), to the 
systems of anionic lipids mixed either with a cationic lipid or a neutral lipid. We have solved 
the model within self consistent field theory, and obtained the phase diagram for these 
systems showing the transition from the lamellar to inverted hexagonal phases which occurs 
at a value of the pH which depends continuously on the membrane composition. Identifying 
the instability of unilamellar vesicles as evidenced by their fusion with the instability of the 
lamellar phase of the lipid mixture, we obtain good agreement between our calculations and 
the experimental data of Hafez et al. for the mixed anionic, cationic system which they 
studied (Hafez et al., 2000). We have interpreted the instability in terms of the counter ions 
which, by increasing the effective size of the anionic lipid, stabilize the lamellar phase of the 
mixture. As cationic lipids are added, counter ions are subtracted. When their number gets 
too low, the vesicle becomes unstable, and the lipids try to revert to an inverted hexagonal 
phase. As noted by Hafez et al., the pH at which this instability occurs is a continuous 
function of the composition of the vesicle, and can therefore be readily tuned to occur at a 
biologically relevant value of pH. Hence such a vesicle has the possibility of being used for 
drug delivery. Applying similar reasoning, we suggested that a vesicle consisting of a neutral, 
non-lamellar forming lipid, like PE, and stabilized by the presence of an anionic lipid could 
also be used to make pH sensitive vesicles whose pH at the point of instability could also 
be readily tuned. We applied our model to such a system and found this to be true, and 
that both sorts of vesicles displayed a region in which the pH at instability was a sensitive 
function of composition, and thus readily tuned. 
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It might be argued that the instabihty of the mixed, charged, hpid vesicles prepared by 
Hafez et ah does not reflect the L^ to Hu phase transition of the bulk system, but rather 
is induced by a phase separation of the components. Such a scenario would be contrary to 
experimental evidence previously cited (Hope et al., 1983; Ellens et al., 1986) and seems to 
us most unlikely. Not only is the mixed system energetically favorable, due to the Coulomb 
attraction of the two components, but it is also entropically favorable because formation 
of the mixture liberates the counter ions needed to neutralize the phase-separated charged 
components (Paulsen et al., 1988). 

Throughout this paper we have reiterated the view that the counter ions, which we have 
treated as a proton with four waters of hydration, H9O4", play an important role in bringing 
about the transition with pH by effectively increasing the volume of the headgroup of the 
anionic lipid. This is reasonable given the well-known stabilization of the lamellar phase with 
respect to the inverted hexagonal one as the headgroup volume increases (Gruner, 1989). We 
further understand the importance of the counter ions as follows. The Coulomb interaction 
in this system has several effects. One of them is that just noted; the charged headgroups 
attract counter ions which have gained an effective volume due to their interaction with the 
dipoles of water, and these tend to stabilize the lamellar phase with respect to the inverted 
hexagonal. However the Coulomb interaction has another effect, which is less appreciated. 
The Coulomb repulsion between the headgroups tends to separate the lipids, just as an 
increase in temperature does, and the consequence is the same; the tails have more room in 
which to move, and this tends to destabilize the lamellar phase with respect to the inverted 
hexagonal one. Thus there is a competition between these two effects, the outcome of which 
depends, among other things, on the effective volume of the stabilizing counter ions. To 
confirm this picture, we have verified by explicit calculation that, if the effective volume of 
the counter ions is too small, the destabilizing tendency of the repulsion between headgroups 
overcomes the stabilizing tendency of the added volume of counter ions, with the result 
that the system makes a transition from lamellar to inverted hexagonal on increasing the 
pH. This is, of course, opposite to experimental observation and to our calculations using 
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counter ions of volume appropriate to those found in water. Nonetheless, the existence of this 
competition between different aspects of the Coulomb interaction emphasizes the importance 
of the counter ions and their volume, and indicates at least one mechanism whereby the pH 
at which the lamellar phase becomes unstable should be expected to depend not only upon 
the system composition, but also upon the species of counter ion and the nature of the 
solvent, as is observed (Seddon, 1990). 

We are grateful to Dr. I.M. Hafez for providing a preprint of his work with S. Ansell and 
P.R. Cullis, and for informative correspondence. This work was supported in part by the 
National Science Foundation under grant number DMR9876864. 
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FIGURE LEGENDS 

Figure 1 Phase diagram in the temperature, T*, pH plane for a system of a single anionic 
lipid, solvent, and counter ions. The pK is that of the anionic headgroup. The volume 
of the headgroup relative to that of the entire lipid is 0.286, similar to that of DOPE, 
the relative volume of the solvent is close to that of water, and the relative volume of 
the counter ions is that of HgOf. 

Figure 2 Phase diagram, at fixed temperature, T* = 0.079 and solvent activity, Zg = 3.425, 
of a mixture of ionizable, anionic lipid, and fully ionized cationic lipid. The volume of 
the headgroup of each lipid relative to that of the entire lipid is the same as that in 
Fig. 1. The relative concentration of the ionizable lipid is denoted G. The solid lines 
show the coexistence obtained from our calculation, the solid circles show the data 
from Hafez et al., 2000, and the dashed line their criterion of vanishing counter ion 
density, Eq. |33| . 

Figure 3 For the system of Fig. 2, the spatial distributions of the various components and of 
the charges is shown for the lamellar phase at coexistence with the inverted hexagonal 
phase near pH— pK=0. Part (a); shown are the volume fractions of the headgroup of 
the anionic lipid, 0/^ , of the tails of the anionic lipid, 'yt(pi , of the headgroup of the 

(2) (2) 

cationic lipid, (f)\ , and of the tails of the cationic lipid, 7t0( , of the solvent, 7^0^, 
and of the counter ions, 7c0c- Part (b); shown are the charge distributions, in units of 
evj^^, of the anionic lipid, pj^ , of the cationic lipid, pj^ , of the counter ions, 0c, and 
the total charge distribution. The points x = and x = D correspond to the centers 
of adjacent solvent regions, with D = 2.94Rg being the lamellar period. 

Figure 4 The same quantities which are shown in Fig. 3a and 3b for the lamellar phase are 
shown here for the inverted hexagonal phase with which the lamellar phase coexists. 
The points x = and x = D correspond to the centers of adjacent tubes, with 
D = 3.03Rg. 
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Figure 5 Phase diagram, at fixed temperature, T* = 0.079 and solvent activity, Zg = 3.425, 
of a mixture of ionizable, anionic lipid, and neutral lipid. The volume of the headgroup 
of each lipid relative to that of the entire lipid is the same as that in Fig. 1. The relative 
concentration of the ionizable lipid is denoted O. The solid lines show the coexistence 
obtained from our calculation. For comparison, the dashed dotted lines show the 
calculated coexistence for the ionizable anionic and fully ionized cation lipid system 
shown previously in Fig. 2 as solid lines. 
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Figure 2: X.-J. Li and M. Schick, Biophys. J., top f|~ 



25 



0.7 



g ^-^ yd) / 

S 0.4- Y^i 



/ 






\ 
\ 



(a). 




Figure 3: X.-J. Li and M. Schick, Biophys. J., top f|~ 



26 






o 
> 



0.8 
0.6 
0.4 
0.2 




- \ 



■■^yA 






,(1) / ^t^t N^ 



(a)J 
/ -\ 




0.2 0.4 0.6 0. 




Figure 4: X.-J. Li and M. Schick, Biophys. J., top f|~ 



27 



0.6- 




Figure 5: X.-J. Li and M. Schick, Biophys. J., top f|~ 



28 



